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Abstract 

A new international Pilot Project for the re- 
determination of the ICRF was initiated by the In- 
ternational VLBI Service for Geodesy and Astrom- 
etry (IVS) in January 2005. The purpose of this 
project is to compare the individual CRF solutions 
and to analyze their systematic and random errors 
with focus on the selection of the optimal strategy 
for the next ICRF realization. Eight radio source 
catalogues provided by the IVS Analysis Centers 
GA, SHAO, DGFI, GIUB-BKG, JPL, MAO NANU, 
GSFC, USNO were analyzed. In present study, four 
analytical models were used to investigate the sys- 
tematic differences between solutions: solid rotation, 
rotation and deformation (IERS method) , and expan- 
sion in orthogonal functions: Legendre-Fourier poly- 
nomials and spherical functions. It was found that 
expansions by orthogonal function describe the differ- 
ences between individual catalogues better than the 
two former models. Finally the combined CRF was 
generated. Using the radio source positions from this 
combined catalogue for estimation of EOP has shown 
improvement of the uncertainty of the celestial pole 
offset time series. 

1 Introduction 

Celestial reference system (CRF) as realized by a set 
of coordinates of selected celestial objects is widely 
used for numerous astronomy navigation, time and 
other measurements. The CRF accuracy and sta- 
bility arc all-important for successful solution of all 
these tasks. For millennia, the CRF was based on 
optical observations and star positions. With es- 
tablishment of new observing technique, very long 
baseline intcrfcrometry (VLBI), much more accurate 
CRF realization became available. In 1998, the CRF 
based on the positions of extragalactic radio sources 
has been adopted by IAU (International Astronomy 



Union) as the fundamental celestial reference frame, 
replacing the FK5 optical frame (Arias et al. 1995, 
Ma et al. 1998). 

After publishing of the first VLBI RSC, attempts 
was made to improve the accuracy of radio-band CRF 
by means of constructing of combined catalogues, as 
it was customary for optical astronomy, where funda- 
mental catalogues served as an international standard 
for astrometry and other measurements on the sky. 
Different methods were used to obtain a combined 
RSC, e.g. Walter (1989a,b), Yatskiv & Kur'yanova 
(1990), Kur'yanova & Yatskiv (1993). Also, till 1995, 
IERS (International Earth Rotation Service) used de- 
rived combined RSC for maintenance of the IERS 
Celestial Reference Frame. 

However, starting from 1996, new CRF realization 
was adopted by the IERS, and further approved by 
the IAU in 1998. The first realization of the ICRF 
was based on the refined analysis of VLBI observa- 
tions made at the NASA Goddard Space Flight Cen- 
ter, USA (Ma et al. 1998). All the 608 radio sources 
included in the ICRF was divided into three groups: 
212 defining sources whose coordinates are supposed 
to be kept in the future realizations to maintain the 
ICRF orientation, 294 candidate sources not suffi- 
ciently monitored, and 102 other sources for improv- 
ing the sky coverage. In 1999 and 2004, two ICRF 
extensions ICRF-Ext.l (Ma 2001) and ICRF-Ext.2 
(Fey et al. 2004) have been published. In those ver- 
sions, the positions of 212 defining sources were kept 
the same as obtained in the first ICRF. It should be 
noted that both ICRF extensions were obtained in a 
manner similar to the first realization, i.e. as a result 
of analysis of the VLBI observations at a single analy- 
sis center. The latest ICRF realization, ICRF-Ext.2, 
is hereafter referred to as ICRF. 

In the end of 2004, joint pilot project of the IERS 
and the IVS (International VLBI Service for Geodesy 
and Astrometry, Schlueter et al. 2002) was initiated 
(Ma 2004, Call for Participation). One of the main 



goals of the project was to seek after possible ways to 
improve the existing ICRF. Large experience accrued 
by the optical astrometry during centuries shows that 
combining catalogues of the star positions have bet- 
ter random and systematic accuracy than individual 
catalogues. In particular, the latter can be affected 
by the systematic errors caused by algorithms and 
software used for data processing. Hopefully, a com- 
bining procedure can be used to mitigate influence of 
errors of individual RSC. Another main goal of the 
Pilot Project is to develop new methods of compar- 
ison of RSC adequate to the modern level of their 
precision and accuracy, in other words, their random 
(stochastic) and systematic errors. 
In this paper, the work was made in four steps. 

1. Analysis of the random and systematic errors of 
individual (input) catalogues, and a choice of the 
most adequate method of representation of the 
systematic differences between catalogues. 

2. Determination of the systematic differences be- 
tween the input catalogues and ICRF. 

3. Construction of a combined catalogue in the 
ICRF system (stochastic improvement of the 
ICRF). 

4. Construction of the final combined catalogue 
(systematic improvement of the ICRF). 

First, we searched after optimal method of repre- 
sentation of the systematic errors of the RSC. Then 
we investigated a possibility of improving the ICRF 
by means of combining individual RSC. Four meth- 
ods of representation of the systematic part of differ- 
ences in the RSC have been examined on the basis 
of comparison of the residual differences in the ra- 
dio source coordinates. Eight individual radio source 
catalogues, and ICRF were used in this study. Af- 
ter the most accurate method had been chosen, it 
was used to compute the systematic differences be- 
tween the individual catalogues and ICRF. Finally, 
these differences were used in the procedure of con- 
struction of combined RSC. Thus obtained combined 
catalogue was tested by means of computation of ce- 
lestial pole offset time series with both combined and 
ICRF RSC. Result of this test showed improvement 
of the scatter of the time series when combined RSC 
is used. 

2 Input catalogues 

Input catalogues used in this study were submit- 
ted by eight IVS Analysis Centers: AUS (Geo- 
science Australia), BKG (Bundesamt fiir Kartogra- 
phie und Geodasie, Germany), DGFI (Deutsches 




Figure 1: Sky distribution for 196 defining (top) and 
all 525 sources (bottom) common for the input cata- 
logues. 

Geodatisches Forschungsinstitut, Germany), GSFC 
(NASA Goddard Space Flight Center, USA), JPL 
(Caltech/NASA Jet Propulsion Laboratory, USA), 
MAO (Main Astronomical Observatory of Na- 
tional Academy of Sciences of Ukraine), SHAO 
(Shanghai Astronomical Observatory, China) , USNO 
(U. S. Naval Observatory, USA). Brief description of 
the input catalogues is given in Table [T] 

Usually, investigation of the systematic differences 
between catalogues is made using some set of refer- 
ence sources common for compared catalogues. Com- 
parison of the lists of radio sources included in the 
input catalogues showed that there are 525 sources 
present in all the catalogues, 196 of them belong to 
the list of 212 ICRF defining sources. This should be 
mentioned that we took into account only the sources 
which have at least 15 observations in 2 sessions. 
After such a selection, the total number of sources 
present in all the input catalogues amounts to 968. 
Fig.[T]shows the distribution of common sources over 
the sky. 

Both source lists, as well as other subsets of 525 
sources, can be used as reference for analysis of the 
systematic differences between catalogues. All the 
computations described below were carried out for 
both 196 and 525 sources. At this paper, we present 
only the results obtained with the first list of 196 
sources. Although definite differences in results were 
found, the main conclusions made in this study do 



2 



Table 1: Input catalogues. The last column shows number of the sources in the catalogue and number of 
reference sources used to tie the orientation of the catalogue to ICRF. 



Center Software 



AUS OCCAM 

BKG Calc/Solve 

DGFI OCCAM 

GSFC Calc/Solve 

JPL MODEST 

MAO SteelBreeze 

SHAO Calc/Solve 

USNO Calc/Solve 



Time span 
month/year 
11/1979 - 12/2004 
01/1980 - 01/2005 
01/1980 - 01/2005 
08/1979 - 01/2005 
10/1978 - 01/2005 
10/1980 - 01/2005 
04/1980 - 01/2005 
09/1979 - 01/2005 



Number of 

delays sources 

3208197 737(207) 

4031453 748(212) 

3650771 686(199) 

4574189 954(212) 

3575847 734(2) 

3773765 685(25) 

4431503 813(212) 

4252684 943(207) 



Aa AtS 



BKG DGFI JPL USNO GSFC MAO SHAO ICRT BKG DGFI JPL USNO GSFC MAO SHAO ITRF 




Figure 2: WRMS differences between the input cat- 
alogues and ICRF. Unit: /ias. 

not depend on the reference source list. 

Weighted root-mean-square (WRMS) differences of 
the radio source coordinates between the input cata- 
logues and ICRF are shown in Fig. [2] (/ias stands for 
microarcseconds) . One can see from Fig. |2j WRMS 
differences have the least values for catalogues com- 
puted with Calc/Solve software, both for intercom- 
parison of these catalogues and their comparison with 
ICRF. The latter most probably is caused by the fact 
that the ICRF was constructed using Calc/Solve soft- 
ware. Large WRMS differences between JPL and 
other catalogues may be caused by its orientation 
to ICRF which has been defined by only two refer- 
ence sources, unlike other catalogues, for which much 
longer lists of reference sources were used. Catalogue 
AUS shows the greatest differences with other cata- 
logues, probably, because it is the only catalogue con- 
structed using the Least Squares Collocation method, 
while other Analysis Centers used conventional Least 
Squares. One can see that the DGFI catalogue, 
also constructed using OCCAM software, but Least 
Squares version, does not stand out against other cat- 
alogues. 

The WRMS differences between the input cata- 
logues themselves may be not valuable, especially 



when they are caused by the difference in orientation 
of the catalogues' axes, which can be easily accounted 
for during combination. More interesting fact is that 
all the input catalogues demonstrate rather large dif- 
ferences with the ICRF, which may indicate signifi- 
cant systematic errors in the latter. This preliminary 
conclusion will be confirmed by the analysis of the 
systematic differences given below. 

3 Analytical representation of 
the systematic differences 

Suppose, we have two catalogues of position of ce- 
lestial sources given as set of spherical coordinates 
ctii,6ii and oi2i-,&ii-, where i is the source number. 
Denote the differences between spherical coordinates 
of the i-th object given in two catalogues as Acti = 
ol\i — ct2i, ASi — &\i — The aim of an analytical 
representation of the systematic differences between 
two catalogues is to obtain formulae 

Aa = f a (a, 8) , 
A5 = f s (a,5), 

(1) 

{ AS } = 

which provides minimum residuals between Aaj, Adi 
observed and computed analytically. In fact, such a 
representation of the systematic differences is a kind 
of a low-pass filter, which allow us to smooth stochas- 
tic errors in coordinates. Having such a representa- 
tion, one can reduce a RSC to the system of another 
catalogue. 

In this section, we compare four methods of an- 
alytical representation of the systematic differences 
between radio source positions given in various cata- 
logues. Those methods are: simple rotation around 
three Cartesian axes hereafter referenced to as "R" , 
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rotation plus deformation used by IERS ("RD"), 
Brosche's method ("B") and expansion in Legcndrc- 
Fourier functions ("LF"). 

It should be noted that in usual astrometric prac- 
tice differences between source positions in right as- 
cension are used as Aacos<5, which reflects the ge- 
ometry of the celestial sphere. However, hereafter we 
use Aa because IERS's method is formulated only for 
this type of differences. Other method listed below 
can be easily adapted to Aa cos 6. 

3.1 Simple rotation 

Two given catalogues realize two Cartesian coordi- 
nate systems X{Y\Z\ and X^Y^Z^- Then differences 
between source positions in two catalogues can be 
represented as the result of rotation of the second co- 
ordinate system with respect to the first coordinate 
system about axes XY Z by three angles A\, A2,A^. 
Then the systematic differences between two cata- 
logues can be expressed as (Walter & Sovers, 2000) 



Aa = A\ tan S cos a + A2 tan 8 sin a — A3 
A5 = —Ai sin a + A2 cos a . 



3.2 Rotation with deformation 



(2) 



This method of representation of the systematic dif- 
ferences between radio source catalogues was pro- 
posed by Arias and Bouquillon (2004), and it is used 
by the IERS since 1995. The authors added to (|) 
three supplement terms to account for some specific 
errors of VLBI catalogues. In this method, the sys- 
tematic differences between two catalogues are ap- 
proximated by 



Aa = Ai tan S cos a + A2 tan 6 sin a — A3 
AS = —A\ sin a + A2 cos a + Dg + Bg , 



3.3 Brosche's method 



D a S, 
(3) 



In this method, as well as in the next one, an ex- 
pansion of the differences in source positions between 
catalogues in orthogonal functions is used. Large ex- 
perience collected by the optical astrometry proved 
that such an expansion provides the highest accuracy 
of the representation of the systematic errors of cat- 
alogues of celestial source positions. In this case, the 
general representation of the differences is given 

by 



where bj are the coefficient to be found from analysis 
of the differences. According to Brosche (1966) 

r p n0 (s), k = 0,1^1, 

Yj(a,S) = i P nk (5)sm(ka), fc^0,Z = 0, (5) 
[ P nk (5)cos(ka), k^ 0,1=1, 

where P n k, associated Legendre polynomials are 
given by 



Pnk(S) = cos k S 



[p/2] {-iy n (p-f) 
i p 5 + -a — — sinP ~ 2AI 5 

m=i n 2^(2n- 2v + 1) 



where p = n — k, [p/2] is entier of p/2, n — 0, 1, 
k = 0, 1, 2, . . . , n, j = n 2 + 2k + I - 1. 



(G) 



3.4 Legendre-Fourier functions 

Bien et al. (1978) has proposed to use another set 
of orthogonal functions for better representation of 
systematic differences between catalogues, especially 
in the polar regions. In this case Q is given by 



Aa 
AS 



b n kiY nk i(a, S) 

nkl 



Y nk i(a,5) = R nkl L n {s\n5)F k i{a) 



Here we omit the Hermite function included in the 
original expression of Bien et al. (1978) to account 
for the source brightness. Commonly speaking, such 
a dependence may do exist in the case of VLBI ob- 
servations too, and it worth separate investigating. 
Legendre polynomials can be computed using recur- 
sion 

L\ = sin 5 . 

L n+ i(sin<5) = 2 r [ l ^i sin 5 L n (sin S) - ^ T i n _i(sin S) , 



(8) 



ra = 2,3,... . 
Fourier functions are given by 

!1, fc = 0, l = -l, 

cos(kla), k^O, 1 = 1, (9) 
sin(— kla), k ^ 0, I = —1 . 

Lastly, normalizing functions are given by 

1, fe = 0, 



R, 



kl 



V2n + 1 



V2, k^0. 



(10) 
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All the four methods described above were applied 
to the differences between each of eight input cat- 
alogues and ICRF for 196 common defining sources 
(see section refinput. For this purpose, the coeffi- 
cients of pj, (pj, ^ were found by means of 
Least Squares adjustment. After that we computed 
the residuals between original differences and those 



((5}, Q. The re- 
and Table [2] Repre- 



computed by formulae (|2]), 
suits are presented in Fig. [3 
scntation with the Brosche's model is not shown in 
Fig. [3] because corresponding surface practically co- 
incides with the Legendre-Fourier expansion. One 
can see that expansion in Legendre-Fourier functions 
provides the least residuals, it i.e. most accurate rep- 
resentation of the systematic differences between cat- 
alogues. Expansion in spherical functions (Brosche's 
method) gives worse accuracy. As to the first two 
methods, they seem to be not adequate to actual er- 
rors of modern RSC. 



4 Combined catalogue in the 
ICRF system 

At the next step, the systematic differences between 
the input catalogues and ICRF found by the LF 
method as described in the previous section were ap- 
plied to all of the input catalogues with the view to 
transforming them to the ICRF system. After that, 
the coordinates of all the sources in transformed cat- 
alogues were averaged with weights depending on the 
formal errors of coordinates. In result, the combined 
catalogue RSC(PUL)07C01 was constructed. This 
catalogue containing all the 968 sources present in 
the input catalogues can be considered as a stochastic 
improvement of the ICRF. Fig. [4] shows the WRMS 
differences between the input catalogues transformed 
to the ICRF system and RSC(PUL)07C01. 

Comparing Fig.|2]and|4] one can conclude that the 
differences between the input catalogues contain not 
only systematic part described by analytical repre- 
sentation, but also, in some cases, significant stochas- 
tic and/or high-frequency components. 

Fig. [5] shows the systematic difference between 
combined catalogue RSC(PUL)07C01 and ICRF. 
One can see that the catalogue RSC(PUL)07C01 rep- 
resents the ICRF system at a level of about 10 ^as. 

5 Final combined catalogue 

Final combined catalogue was constructed in the fol- 
lowing way. Let us call the system of a catalogue the 
set of coefficients of Legendre-Fourier functions ob- 
tained for given catalogue. Thus we have eight sys- 













Figure 3: Analytical representation of the differences 
between the input catalogues and ICRF: R (dark 
grey), RD (grey), LF (light grey). Original differ- 
ences are shown in black lines. Horizontal axes show 
right ascension (right) and declination (left) in de- 
grees. Unit: (ias. 
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Table 2: WRMS residuals between the input catalogues and ICRF before (Raw) and after approximation of 
the systematic differences (see notation of the methods in text). Results related to the LF method providing 
the best approximation arc shown in bold. Unit: fia,s. 





AUS 


BKG 


DGFI 


JPL 


USNO 


GSFC 


MAO 


SHAO 












Aa 








Rciw 


304 


277 


343 


331 


253 


277 


294 


273 


R 


301 


271 


342 


308 


249 


274 


286 


271 


RD 


299 


265 


342 


308 


247 


273 


285 


270 


B 


170 


177 


237 


238 


172 


191 


203 


193 


LF 


106 


125 


164 


172 


122 

AS 


144 


152 


145 


Raw 


337 


225 


284 


277 


221 


238 


251 


240 


R 


337 


225 


284 


273 


221 


238 


251 


240 


RD 


333 


224 


283 


273 


221 


237 


251 


239 


B 


180 


159 


178 


182 


152 


158 


169 


166 


LF 


111 


109 


112 


127 


104 


106 


134 


111 



Act AS 



BKG DGFI JPLLSNO GSFC MAD SHAO H'RFCOMB BKG DGFI JPL USNO GSFC MAO SHAO ICRF COMB 




Figure 4: WRMS differences between the input cata- 
logues transformed to the ICRF system and combined 
catalogue RSC(PUL)07C01. Unit: ^as. 



Ac* AS 




Figure 5: Differences between RSC(PUL)07C01 and 
ICRF. Horizontal axes show right ascension (right) 
and declination (left) in degrees. Unit: /xas. 




Figure 6: WRMS differences between the input cat- 
alogues and combined catalogue RSC(PUL)07C02. 
Unit: /ias. 

terns corresponding to the eight input catalogues. At 
the first iteration, these systems were averaged with- 
out weights. Then the WRMS differences between 
this average system and the systems of the input cat- 
alogues for the bins 10° (a) x 5° (S) were computed. 
Thus obtained WRMS were used for weighting of the 
input catalogues at the second iteration. Although 
we use different weights of the input catalogues for 
different bins, they are close for each catalogue. Fi- 
nal weights of the catalogues averaged over the sky 
are given in Table [3j 

Thus computed average system was added to 
the first combined catalogue RSC(PUL)07C01. In 
result, final combined catalogue, RSC(PUL)07C02 
have been obtained. It can be considered as both 
stochastic and systematic improvement of the ICRF. 
In Fig. [6] WRMS differences between the input cata- 
logues and combined catalogue RSC(PUL)07C02 are 
shown, and Fig. [7] presents the systematic differences 
between the input catalogues and RSC(PUL)07C02. 

Comparison of RSC(PUL)07C02 and ICRF is pre- 
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Table 3: Weights of the input catalogues applied during combination, averaged over the sky. 



AUS 


BKG 


DGFI 


GSFC 


JPL 


MAO 


SHAO 


USNO 


a 0.246 
6 0.205 


0.671 
1.220 


0.464 
0.586 


1.993 
1.921 


0.254 
0.446 


0.558 
0.541 


1.792 
1.459 


2.062 
1.927 



Table 4: WRMS differences between celestial pole 
offset series computed with two CRF realization and 
IAU2000A model corrected for FCN contribution. 
Unit: fias. 



Catalogue 


x c 


Y c 


ICRF-Ext.2 


103 


100 


RSC(PUL)07C02 


98 


98 



sented in Fig.[8j which show the result of expansion of 
the differences between two catalogues in Legendre- 
Fouricr functions. Results of this comparison lead us 
to the supposition that ICRF may have significant 
systematic errors. 

6 Comparison with observa- 
tions 

It is important to assess an actual accuracy of ob- 
tained catalogue (as well as other CRF realizations). 
Unfortunately, existing methods of comparison of 
catalogues allow us to investigate only differences be- 
tween catalogues. Here, we use a test, which can help 
us to get some independent estimate of the quality of 
our combined catalogue. For this purpose, we com- 
pute two celestial pole offset time series from process- 
ing of Rl and R4 IVS observing programs observed 
in the period 2002-2006 with two radio source cata- 
logues, ICRF-Ext.2 and RSC(PUL)07C02. Then we 
compute the WRMS differences between computed 
celestial pole offsets and Free Core Nutation model. 
The result of this test presented in Table [4] shows 
clear improvement of the scatter of celestial pole off- 
set estimates. 

7 Conclusion 

In this paper, we have constructed a new combined 
catalogue of radio source coordinates. For this study 
we used eight catalogues submitted by IVS Analy- 
sis Centers in the framework of the IERS/IVS Pilot 
Project on the future realization of ICRF. 

First, we have examined four methods of analyti- 
cal representation of systematic differences between 
catalogues of radio source coordinates. Representa- 



tion by means of expansion in Legendre- Fourier func- 
tions is proved to be the most accurate method. Two 
methods usually used for comparison of radio source 
catalogues, simple rotation and rotation with defor- 
mation seem to be not suitable for investigation of 
modern radio source catalogues, and should be re- 
placed by more adequate one. 

Two combined radio source catalogues have been 
constructed. The first of them is obtained as weighted 
average of the input catalogues corrected for sys- 
tematic differences with ICRF. It can be considered 
as stochastic improvement of the current realization 
of ICRF. Second combined catalogue have been ob- 
tained from the first one after applying the weighted 
average systematic difference between the input cat- 
alogues and ICRF, which allows one to account also 
for possible ICRF systematic errors. 

To compare our combined catalogue with ICRF, we 
used two tests, which allow us to independently es- 
timate the scatter of celestial pole offsets time series 
obtained from processing of VLBI observations us- 
ing ICRF and combined catalogue. Both tests have 
showed improvement of the results. 

The results obtained in this paper allow us to make 
a conclusion that ICRF-Ext.2 may have significant 
systematic errors, most probably caused by fixing the 
coordinates of 212 defining sources to its values ob- 
tained in the first ICRF version of 1995. 

Further development of this study may include: 

• More detailed analysis of stochastic and system- 
atic errors of radio source catalogues. 

• Estimation of possible impact of the high- 
frequency systematic errors in source position on 
the orientation of the catalogue axes. 

• Analysis of the reasons of systematic differences 
between radio source catalogues. 

• Careful selection of the input catalogues and ref- 
erence sources. 

• Elaboration of weighting method. 

• Development of more sensitive methods of the 
assessment of the accuracy of CRF realizations. 

It is also interesting to compare different methods 
of construction of combined catalogue of radio source 
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Figure 7: Systematic differences between the 
input catalogues and combined catalogue 
RSC(PUL)07C02. Horizontal axes show right 
ascension (right) and declination (left) in degrees. 
Unit: ^,as. 



Aa AS 




Figure 8: Differences between RSC(PUL)07C02 and 
fCRF. Horizontal axes show right ascension (right) 
and declination (left) in degrees. Unit: fias. 

coordinates, in particular a method used so-called arc 
approach developed by Yatskiv & Kur'yanova (1990). 
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